RFoptions(seed=0) ## *ANY* simulation will have the random seed 0; set
## RFoptions(seed=NA) to make them all random again
z <- 1:4
x <- cbind(z, 0)
y <- cbind(0, z)
model <- RMwhittle(nu=0.5)
RFcov(model, x, y, grid=FALSE)## standard is (cartesian) models
## same as above, but explicite:
RFcov(model, x, y, grid=FALSE, coord_sys="cartesian")
## model is valid not on a sphere; x,y coordinates are
## transformed from earth coordinates to sphereical coordinates
RFcov(model, x, y, grid=FALSE, coord_sys="earth")
## now comparable the scale chosen sucht that the covariance
## values are comparable to those int the cartesian case
RFcov(RMS(model, s= 1 / 180 * pi), x, y, grid=FALSE,
coord_sys="earth")
## projection onto a plane first. Then the scale is interpreted
## in the usual, i.e. cartesian, sense:
RFoptions(zenit = c(2.5, 2.5))
RFcov(model, x, y, grid=FALSE,
coord_sys="earth", new_coord_sys="orthographic")
## again, here the scale is chosen to comparable to cartesian case
## here the (standard) units are [km]
RFcov(RMS(model, s= 6350 / 180 * pi), x, y, grid=FALSE,
coord_sys="earth", new_coord_sys="orthographic")
## as above, but in miles
RFcov(RMS(model, s= 3750 / 180 * pi), x, y, grid=FALSE,
coord_sys="earth", new_coord_sys="orthographic",
new_coordunits="miles")
Run the code above in your browser using DataLab